0.1 Direct observation of the effect of isotope-induced- 
disorder on exciton binding energy in LiH^Di ^ mixed 

. . . crystals. 

00 

o 

O ■ V.G. Plekhanov 

(N 

Computer Science College, Erika 7a, Tallinn, 10416, 
ESTONIA 

Abstract. The results of the quantitative study of the renormalization of 

binding energy the Wannier - Mott exciton by the isotope effect are present 

for the first time. For this purpose accurate measurements of the intrinsic 

O . luminescence and mirror reflection spectra of LiH^^Di-j: mixed crystals with 

a clean surface in the temperature range 2 - 100 K were carried. Nonlinear 
dependence of exciton binding energy E;, on the isotope mass Eb ^ f(x) is 
caused by the bands fluctuation broadening which is connected with the isotope- 
induced-disorder. Temperature dependence of exciton binding energy is briefly 
discussed. The extrapolation of the asymptotic linear behavior of the exciton 

a maximum energy to T = K enables to estimate the zero - point renormalization 

of the exciton binding energy. 
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A wide variety of novel isotope effects have been discovered on last four 
^ I decades [1-6] owing to the availability of high-quality bulk semiconductor and 

insulator crystals with controlled isotopic composition (see, also, reviews [7-10]). 
Recent high resolution spectroscopic studies of excitonic and impurity transi- 
ts] ' tions in high-quality samples of isotopically enriched Si have discovered the 

broadening bound excitons emission (absorption) lines connected with isotope- 
f~^ ' induced-disorder as well as the dependence of their binding energy on the isotope 

00 ! mass [11-13]. The last effect was early observed on the bound excitons in dia- 

^D • mond [14-15] and more earlier on the free excitons in LiH2:Di_2; mixed crystals 

[16]. 

As is well-known (see, for example [16, 17, 4]) the band gap energy Eg in 
the T — >0 limit has a dependence on the average isotopic mass M due to two 
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j^ ' effects: a) the renormalization of Eg by the electron-phonon interaction coupled 

with the dependence of the zero-point amplitudes on M (see, also [18]) and b) 
the dependence of the lattice constant on M, leading to a change in Eg through 
the hydrostatic deformation potential. The electron-phonon term is dominant 
[19, 20] and in the case of semiconductor crystals (C; Ge; Si) with a weak isotope 

1/2 

scattering potential is varied approximately as M . The value of the T =0 
electron-phonon renormalization energy contribution to Eg can be independly 
determined from an extrapolation of the high temperature linear dependence 



of Eg ~ f(T) to T=0 but again it should be remembered that Eg ^ f(T) also 
has a small contribution from change in volume with temperature. Therefore 
we should distinguish between effects of the average atomic mass (which imply 
the virtual crystal approximation - VGA) and effects of the mass fluctuations 
(randomness) superimposed onto the virtual crystal. 

In this communication we report the first results of the quantitative study 
the dependence of free exciton's binding energy on isotope mass as well as on the 
temperature in LiHa;Di_a; mixed crystals. We'll show that the exciton binding 
energy increase by 10±1 meV from ^LiH to ^LiD. Moreover, the dependence of 
Efc ~ f(x) is nonlinear. The last effect, as will be shown below, is caused by 
the isotope-induced-disorder of LiH2;Di_2; mixed crystals. Briefly part of these 
results have been published in [21]. 

Specimens of LiH, LiD LiHa;Di_a; (0< x < 1) as well as LiHa;Fi_a; were grown 
from melt using the modified Bridgman - Stockbarger method (see, also [22]). 
This technique is described many times early (see, e.g. [18, 2]). To improve the 
stoichiometric composition with respect to hydrogen (deuterium), the crystals 
grown were additionally annealed in an atmosphere of hydrogen or deuterium 
at a gas pressure of 3 to 5 atm. and a temperature 500 to 550°C (melting point 
is 961 and 964 K for LiH and LiD, respectively). For some crystals the time of 
heat treatment was as long as 20 days (for more details see [2,7]). 

Given the high reactivity of freshly cleaved LiH crystals in the atmosphere, 
we had to develop a procedure of cleaving which would not only ensure an 
initially clean surface, but also allow us to keep it fresh for a few hours (the time 
taken to complete an experiment). These requirements are satisfied by the well- 
known method of cleaving directly in the helium chamber of the optical cryostat 
under liquid or superfluid helium, first tried in Ref. [23] . We did not notice any 
changes in the spectra of refiection or luminescence while working for 10 tO 16 
hours with surface prepared in this way. The device for cleaving the crystals 
had three degrees of freedom and rotated through 90°, which greatly helped in 
carrying out the experiments. As a rule, specimens for experiments were cleaved 
off bulk high - quality crystals. The experimental setup for measuring the low - 
temperature reflection and luminescence spectra has been described more than 
once (see, e.g. [7,18]) and consists of a double grating or prism monochromators, 
an immersion helium cryostat, and a photovoltaic detector (in photon counting 
mode). The results presented in this paper were obtained from a clean crystal 
surface cleaved, as described above, directly in the bath of helium cryostat. 
The spectra of mirror reflection were measured using an the angle of incidence 
of 45*^. For our studies we selected specimens which exhibited low dependence 
of exciton spectra of refiection and luminescence on surface features. 

As demonstrated earlier (see, e.g. [7]) most low-energy electron excitations 
in LiH (LiD) crystals are large-radius excitons. The spectrum of exciton pho- 
toluminescence of LiH crystals cleaved in liquid helium consists of a narrow (in 
the best crystals, its half-width is the AE ^ lOmeV) phononless emission line 
and its broader phonon repetitions, which arise due to radiative annihilation of 
excitons with the production of one to five longitudinal (LO) phonons (see. Fig. 
1). The phononless emission line coincides in an almost resonant way with the 



reflection line (see, insert Fig. 1) of the exciton ground state which is another 
indication of direct electron transition (Xi - X4 [2]). The lines of phonon repli- 
cas form an equidistant series biased toward lower energies from the resonance 
emission line of excitons. The energy difference between these lines is about 
140 meV, which is close to the calculated energy of the LO phonon in the mid- 
dle of the Brillouin zone [21] and measured in [19]. The isotopic shift of the 
zero-phonon emission line of LiH crystals equals that in reflection spectra, 103 
meV (see, also, insert in Fig.l). As the deuterium concentrations increases, the 
width of the long-wave maximum in reflection (and the width of the phonon- 
less line in luminescence) spectra broadens and maximum shifts towards the 
shorter wavelengths. As can clearly see in Fig. 1, all spectra exhibit a similar 
long-wave structure. This circumstance allows us to attribute this structure to 
the excitation of the ground (Is) and the first excited (2s) exciton states [2]. 
Three effects are distinctly shown in the reflection spectra by an increase in 
deuterium concentration: 

1. The shortwavelength shift of the reflection (as luminescence) spectrum as 
a whole; 

2. The different velocity shift of the exciton maximum of the ground and 
the first excited states on the temperature (see. Fig. 2); 

3. Broadening of the long- wavelength maximum due to excitation of the 
ground exciton state. 

By the way we should note that the first excited exciton state is very clearly 
seen in the luminescence spectra too (see Fig. 3). Although two distinct contri- 
bution to AEg (and respectively Eb) are present that due to the effect of the 
vibration on the lattice constant and the direct effect of the electron-phonon 
interaction [19, 20] for the present task of the dependence of the exciton bind- 
ing energy on the isotope mass it will be considered only the mechanism of 
exciton-phonon interaction. 

In accordance with the second effect in reflection spectra (see above point 

2), as is to be expected, there is nonlinear dependence of the exciton binding 

energy (E^ = 4/3A12) on the isotope concentration. Displayed curve in Fig. 

4 is the Wannier - Mott exciton binding energy value which are determined 

4 
by the hydrogenlike expression E(, ~ ^^^ (^^^ using for this Ei and E2 values 

from the reflection and luminescence spectra (see Figs. 1 and 3). The nonlinear 
nature of this dependence is similar to the theoretical results [25] , where the 
influence of the chemical disorder of the crystal lattice on the Wannier - Mott 
exciton binding energy was seen and obtained a qualitative (not quantitative!) 
agreement with experimental results on the GaAsa^Pi-^, [26] mixed crystals. 

Before the comparison of our experimental results with the theory developed 
by Elliott and Kanehisa [25], it would be prudent to briefly review main prop- 
erties of their theoretical model. One of the principal result of paper [25] is the 
nonlinear dependence of exciton binding energy E^ on the concentration. As a 
consequence, the binding energy at half-and-half concentrations is less than the 
value derived from the crystal virtual model. According to Ref. [25] this model 
considers an exciton with a direct gap of a semiconductor alloy. Such a system 



consists of an electron (particle 1) in the conduction band (c) with mass nic and 
a hole (particle 2) in the valence band (v) with mass mv The problem of the 
exciton in disordered systems is to solve the Hamiltonian 

H = ^V2nic +^V2mv+ uCti - V2) + Vc(Yi) + ¥,(^2), (1) 
with both the Coulomb interaction u and the potential V^ due to disorder 
(v — c,v). Reference [25] neglected disorder - induced interband mixing. As 
it is well known, in place of the electron- hole coordinates, ( r 1, p 1) and ( r 2, 
P2), one may introduce the center-of-mass and relative coordinates, (R, P) 
and (V, 1?) to rewrite (1) as 

H = 1?V2m, + u(y) + 1^V2M + Vc(R + mvY/M) + Vv(R - mcV/M), 

(2) 

where /i^ and M are the reduced and total masses of cxcitons, respectively. 
Because of the random potential, the translational and relative degrees of free- 
dom cannot be decoupled. This is essentially difficult when considering the two- 
body problem in a disordered system (see [25] and references therein). However, 
when the exciton state in question is well separated from other states so that 
the energy spacing is much larger than the translational width and disorder, 

one can forget about the relative motion (H^ = ^^ -I- u( r )) and just apply 
any single-particle alloy theory solely to their translational motion. For each 
exciton state the translational part of Hamiltonian in this case is 

Hi = _P V2M_+ Vc(R) + Vv(R) (3). 

Here Vc and Vv are averages of Vc and Vv with respect to the relative state 
(/), for example: 

K{R) = Jd^V 10 (V)|2Ve[i? + ^V] (4). 

This approach is very similar to the Born-Oppenheimer adiabatic approxi- 
mation. Such situations hold in some mixed alkali halide crystals and probably 
A2B6 crystals. On the contrary, when the exciton binding energy is compa- 
rable to the disorder energy, the adiabatic approximation breaks down, and 
it is essential to take into account the effect of disorder on both the transla- 
tional and relative motions. This is the case with the Wannier-Mott exciton 
in A3B5 alloys, for which the Elliott and Kanehisa model was developed. In 
this case the solution task is to start from the independent electron and hole by 
neglecting u in (2) and then to take into consideration the Coulomb interaction 
between the average electron and average hole. A further simplified approach 
adopted in the literature (see, for example [27] and references therein) in solving 
the Bethe-Salpeter equation [28] is to suppose a free-electron-like one particle 
Green's function with a built-in width to allow for the random potential due to 
disorder. In the cited theoretical model [25], the average (or "virtual crystal") 
gap is given by 

E-(x) = Eo + ((5c-<^v)(x-l/2), (5) 

where Eo is average gap, (5c, S^ are the values of the fluctuation broadening of 
the conduction and valence bands, respectively. Reference [25] also assumed the 
Hubbard density of states for both the conduction and valence bands with width 
Wc and Wv, respectively, as well as similar dispersion in both bands. With this 



assumption the exciton binding energy has been calculated according to the 
coherent potential approximation CPA model. As is well-known, the main idea 
of the coherent potential methods in the introduction of an auxiliary medium 
with a regular , i.e., spatially periodic potential. By definition in the model 
this potential is also complex. The formalism of the coherent potential model, 
convenient for performing calculation and does not include fitting parameters, 
because of using of the density of phonon states from virtual crystal model 
which used virtual crystal approximation (VGA). It should be added here the 
key feature of the model developed in Ref. 25 is the short-range nature of the 
Coulomb potential (for details, see e.g. [7, 29]). 

The data from Fig. 1 and other published sources [2, 7, 21] were used for 
plotting the energy Eb as a function of isotopic concentration x in Fig. 4. The 
binding energy (defined as the band edge minus the exciton energy) is given by 
[25]: 

K'' = Uo + 2|r + W (6). ^ 

In the last relation Uo is the coupling constant at the total exciton momen- 
tum g = 0. 

Theoretical description of the binding energy of Wannicr- Mott excitons as 
a function of concentration x was based on the polynomial derived by Elliott 
and cowokers [25]: 



Eb - E"^^- Eb 
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Eeff, (7) 



Eeff = X (1-x) ^, (8) 

where W ~ Wc + Wv, and Wc and Wv are the widths of the conduction 
band and the valence band which are equal to 21 eV [30] and 6 eV [31, 32] 
respectively. Here Ebow is the curvature parameter found from the function Eg 
oc f(x) {Fjbow — 0.046 eV [7]); 5c and 5^ are the magnitudes of the fluctuation 
smearing of the valence band and the conduction band edges , 5c = 0.103 eV 
and 5v = - 0.331 eV. As follows from Fig. 4, these values of the parameters give 
a good enough description of the nonlinear dependence of the binding energy of 
Wannier-Mott exciton in disordered medium isotope - mixed crystals LiH^iDi-a;. 
This agreement between theory and experiment once again proves the inherent 
consistency of the model proposed by Kanehisa and Elliott, since the isotopic 
substitution affects the short-range part of the interaction potential. 

In this way, the nonlinear dependence of the binding energy of Wannier-Mott 
exciton is caused by isotopic disordering of the crystal lattice. As is seen from 
Fig. 4 the exciton binding energy decreasing (relative linear law (VGA) - see 
dashed line in Fig. 4) in the vicinity of the middle meaning concentration really 
calls out the fluctuated broadening of the edge of the conduction and valence 
bands. In accordance with the theoretical model the last reason gives rise to the 
reduced Eg and there by the shallowing of the exciton levels and, respectively, 
the reduction of Eb. 

As follows from Fig. 1, the addition of deuterium leads not only to the 
short-wave shift of the entire exciton structure (with different rates for Is and 
2s states), but also to a significant broadening of the long- wave exciton reflec- 
tion line. This line is broadened 1.5 - 3-fold upon transition from pure LiH to 



pure LiD. The measure of broadening was the halfwidth of the hne measured in 
the standard way (see e.g. [33]) as the distance between the maximum and the 
minimum in the dispersion gap of the reflection spectrum, taken at half-height. 
The concentration dependence of the halfwidth (AE^) of the long-wave band 
in the exciton reflection spectrum at 2 K is shown in Fig. 5. Despite the large 
spread and the very limited number of concentrations used, one immediately 
recognizes the nonlinear growth of AE^ with decreasing x. A similar concen- 
tration dependence of AE^ in the low-temperature reflection spectra of solid 
solutions of semiconductor compounds A2B6 and A3B5 has been reported more 
than once (see e.g. the review of Elliott and Ipatova [34] and references therein). 
The observed broadening of exciton lines is caused by the interaction of excitons 
with the potential of large-scale fluctuations in the composition of the solid so- 
lution. Efros and colleagues (see e.g. [35]) used the Lifshitz method of optimal 
fluctuation [36] to express the formula for the concentration dependence of the 

broadening of exciton reflection lines: 
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AE^ = 0.5a 






(9) 

where N - the concentration of sublattices nodes where the isotope substi- 
tutes are placed, a — dEg/ dx; rex is the exciton radius which varies from 47 A 
to 42 A upon transition from LiH to LiD [2] . The results of calculation according 
to Eq. (9) are shown in Fig. 5 by a full curve. 

The experimental results lie much closer to this curve than to the straight 
line plotted from the virtual crystal model. At the same time it is clear that 
there is only qualitative agreement between theory and experiment at x > 0.5. 
Nevertheless, even this qualitative analysis clearly points to the nonlinear de- 
pendence of broadening on the concentration of isotopes, and hence to the iso- 
topic disordering. Since isotopic substitution only affects the energy of optical 
phonon for the first time, and, as a consequence, the constant of exciton-phonon 
interaction (in the first place, the Frohlich interaction gp), the nonlinearity of 
functions AEb oc f (x) , AE'^ oc f(x) is mainly related to the nonlinear behavior of 
g| oc f(x) . In this way, the experimental study of the concentration dependence 
of the exciton-phonon interaction constant may throw light on the nature and 
mechanism of the large-scale fluctuations of electron potential in isotopically 
disordered crystals. 

Returning to the results of Fig. 2, let us add that the different temperature 
dependence of exciton peaks of n = Is and 2s exciton states leads to the 
temperature dependence of the binding energies of Wannier - Mott excitons 

Eb - f (T) (10). 

This problem has not received any adequate treatment. More specifically, 
the energy of the exciton binding Ef, in LiH crystals (as well as in mixed crys- 
tals LiHa,Fi_j. (LiDj.Fi_x )) decreases with increasing temperature, whereas Ef, 
increases for excitons of the green and yellow series in CU2O crystals [34]. A 
linear approximation of the exciton binding energy in LiD0.995F0.005 (see curve 
3, Fig. 2) representing Eb at T=0 K gives E(,(0) = 55 meV. From this value we 
can see that renormalization of the binding energy by the zero-point vibrations 
equals approximately ~ 10 % from this value, that is, on the other hand, com- 



poses only half of renormalized exciton binding energy by isotope effect (~ 10 
meV, see above) . It is not excluded that the other part of renormalized exciton 
binding energy is caused by exciton-polar phonon interaction. 

In conclusion, the exciton luminescence and reflection spectra are used in a 
quantitative study of the isotopic and temperature effects in LiHj^Di-^, mixed 
crystals with a clean surface. It was shown that the short-range character of the 
potential of a disordered crystal lattice with isotope substitution is responsible 
for the broadening of the valence and conduction bands. Nonlinear dependence 
of the exciton binding energy on the isotope mass Ef, ~ f(x) is due the isotope- 
induced-disorder of LiIIj;Di_j: mixed crystals. Temperature dependence of ex- 
citon binding energy is briefly discussed. The extrapolation of the asymptotic 
linear behavior of the exciton maximum energy to T = K enables to estimate 
the zero - point renormalization of the exciton binding energy. 
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Figure captions. 

1. Fig. 1. Luminescence spectra of free excitons at 2 K in LiH and LiD 
crystals cleaved in liquid helium. In insert: mirror reflection spectra of crystals. 
Curve 1: LiH; curve 2: LiH2:Di_2: and curve 3: LiD. Curve 4 is the light source 
without crystals. 

2. Fig. 2. Temperature dependence of the distance between the long- 
wavelength peaks (A12) in specular reflection spectra of pure and mixed crystals: 
1 - LiH; 2 - LiD; 3 - LiD0.995F0.005. 



3. Fig. 3. The reflection (1) and luminescence (2) spectra of LiD crystal at 
2 K. 

4. Fig. 4. Concentration dependence of the binding energy of a Wannier 
- Mott exciton at 2 K in LiHj;Di_j; mixed crystals: 1 - VGA approximation 
model; 2 - calculation according to equation (8); experimental points indicated 
by triangles. 

5. Fig. 5. Concentration dependence of the half- width of the ground state 
line of the exciton in the mirror reflection spectrum at 2 K. 1 - VCA approx- 
imation model; 2 - calculation according to equation (9); experimental points 
are indicated by crosses. 



